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B — > iTTr and related decays are studied in the heavy quark limit of QCD using the soft collinear 
effective theory (SCET). We focus on results that follow solely from integrating out the scale mj, 
without expanding the amplitudes for the physics at smaller scales such as ct a (^/ E^Aqcd)- The 
reduction in the number of hadronic parameters in SCET leads to multiple predictions without the 
need of SU(3). We find that the CP-asymmetry in B~ —* ir°K~ should have a similar magnitude 
and the same sign as the well measured asymmetry in B° — > tv + K~ . Our prediction for Hr(K + n~) 
exceeds the current experimental value at the 2a level. We also use our results to determine the 
corrections to the Lipkin and CP-asymmetry sum rules in the standard model and find them to be 
quite small, thus sharpening their utility as a tool to look for new physics. 



I. INTRODUCTION 

Two body nonleptonic decays are the most widely used 
processes to study CP violation in the B system. Due 
to the large mass of the B-meson there is a plethora of 
open channels, each of which provides unique ways for 
testing the consistency of the standard model. For each 
channel obscrvables include the CP averaged branching 
ratios (Br), direct CP asymmetry (Acp = —C), and for 
certain neutral B decays, the time dependent CP asym- 
metry (S). For the decays we are interested in 
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where A is the amplitude of the decay process A — 
A{B — > M1M2), A is the amplitude for CP-conjugate 
process, and q/p is the mixing parameter for B° — B° 
and/or K° — K° mixing. The other parameters in Eq.lfT|) 
are \p\, the final meson momentum in the B rest frame, 
s, a possible identical particle symmetry factor, and Am, 
the difference between mass eigenstates in the neutral B 
two-state system. 

Using the unitarity of the CKM matrix to remove top- 
quark CKM elements, the amplitude for any decay can 
be written with the CKM elements factored out as 
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Theoretical predictions for the 



observables in are often hampered by our ability to 
calculate A UiC . In general the CP-asymmetries depend on 
the ratio of amplitudes |A„/^4 C | and their relative strong 



phase 5. In fact Aqp oc sin(<5), and so non- negligible 
strong dynamics are required for the existence of a direct 
CP asymmetry. 

The parameters A u and A c are in principle different 
for each decay channel. In order to accurately deter- 
mine A u and A c we need model independent methods 
to handle the strong dynamics in these decays. All such 
methods involve systematic expansions of QCD in ratios 
of quark masses and the scale A ~ Aqcd associated with 
hadronization. This includes flavor symmetries for the 
light quarks, SU(2) and SU(3), from m q /A <g; 1, as well 
as expansions for the heavy 6-quark from A/ raj, 1. 
For nonleptonic decays to two light mesons with ener- 
gies E m ~ mb/2, kinematics implies that we must also 
expand in A/ Em « 1. A formalism for systematically 
expanding QCD in this fashion is the soft-collinear ef- 
fective theory (SCET) 0. In nonleptonic B-decays the 
expected accuracy of these expansions are 



SU(2) 
SU(3) 
SCET 



mu,d 

A 
n 
) — 

A 2A 



0.03 < 1 , 



m„ W 771s 

~A ' X ' 



(3) 



E 



m b 



0.2 < 1. 



The flavor symmetries SU(2) and SU(3) provide ampli- 
tude relations between different nonleptonic channels, 
thereby reducing the number of hadronic parameters. 
The expansion in A/m-b ~ A/ Em also reduces the num- 
ber of hadronic parameters. In this case the expansion 
yields factorization theorems for the amplitudes in terms 
of moments of universal hadronic functions. 

In this paper we study standard model predictions for 
B — * Kir, KK, and 7T7t decays. These channels pro- 
vide 25 observables, of which 19 have been measured or 
bounded as summarized in Table [I] We make use of the 
expansions in Eq. ©, focusing on SCET. Our goal is 
to quantify the extent to which the current data agrees 
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A C p = -C 


S 


tt + tt 


5.0 ±0.4 


0.37 ±0.10 


-0.50 ±0.12 


7T%° 


1.45 ±0.29 


0.28 ±0.40 




TT+TT 


5.5 ±0.6 


0.01 ±0.06 


- 


7T-^° 


24.1 ± 1.3 


-0.02 ± 0.04 


- 


7v°K- 


12.1 ±0.8 


0.04 ±0.04 


- 


ir + K~ 


18.9 ±0.7 


-0.115 ±0.018 


- 


7T K 


11.5 ± 1.0 


-0.02 ±0.13 


0.31 ± 0.26 


K+K~ 


0.06 ±0.12 






K°K° 


0.96 ±0.25 






K°K~ 


1.2 ±0.3 







TABLE I: Current B -> tttt, Kir, and KK data 
The S for ttK is S(ir°K s ). 



or disagrees with the standard model in the presence 
of hadronic uncertainties, and to provide a roadmap for 
looking for deviations in future precision measurements 
of these decays. 

The SU(2) isospin symmetry is known to hold to a 
few percent accuracy, and thus almost every analysis 
of nonleptonic decays exploits isospin symmetry. (Elec- 
troweak penguin contributions are simply AI =1/2 and 
AI = 3/2 weak operators, and are not what we mean by 
isospin violation.) Methods for determining or bounding 
a (or 7) using isospin have been discussed in 0,0 and are 
actively used in B — > 7T7t and B —> pp decays. In B — > pp 
this yields a pp = 96° ± 13° Q For B — > 7T7t this analysis 
has significantly larger errors, since the A c amplitudes 
are larger and the asymmetry C(7r°7r ) is not yet mea- 
sured well enough to constrain the hadronic parameters. 
Isospin violating effects have been studied in Q. For 
B — > Kir and B — > KK an SU(2) analysis is not fruit- 
ful since there are more isospin parameters than there are 
measurements, so further information about the hadronic 
parameters is mandatory. 

In B — > 7T7T, even if C(7r°7r°) were known precisely it 
would still be important to have more information about 
the amplitudes A u and A c than isospin provides. For 
example, isospin allows us to test whether 7 7r7r differs 
from the value obtained by global fits [Hj, [13] , 
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However, a deviation in 7 is not the only way that new 
physics can appear in B — > 7T7T decays. Simply fitting 
the full set of SU(2) amplitudes can parameterize away a 
source of new physics. For example, Ref. [13J has argued 
that it is impossible to see new physics in the (7T7r)/ = o 
amplitudes in an isospin based fit. Thus, it is important 
to consider the additional information provided by SU(3) 
or factorization, since this allows us to make additional 
tests of the standard model. The expansion parameters 
here are larger, and so for these analyses it becomes much 
more important to properly assess the theoretical uncer- 
tainties in order to interpret the data. 



The analysis of B — > Kir decays has a rich history 
in the standard model, provoked by the CLEO measure- 
ments [HI that indicated that these decays are domi- 
nated by penguin amplitudes that were larger than ex- 
pected. The dominance by loop effects makes these de- 
cays an ideal place to look for new physics effects. Some 
recent new physics analyses can be found in Refs. |15| . 
This literature is divided on whether or not there are 
hints for new physics in these decays. The main obsta- 
cle is the assessment of the uncertainty of the standard 
model predictions from hadronic interactions. 

Several standard model analyses based on the limit 
m s /A <C 1 fie SU(3) symmetry") have been reported re- 
cently O [H H [H M UK (see also IH H 13 for 
earlier work). In the AS = 1 decays the electroweak 
penguin amplitudes can not be neglected, since they are 
enhanced by CKM factors. Unfortunately the number 
of precise measurements makes it necessary to introduce 
additional "dynamical assumptions" to reduce the num- 
ber of hadronic parameters beyond those in SU(3). In 
some cases efforts are made to estimate a subset of the 
SU(3) violating effects to further reduce the uncertainty. 
The dynamical assumptions rely on additional knowledge 
of the strong matrix elements and in the past were mo- 
tivated by naive factorization or the large N c limit of 
QCD. Our current understanding of the true nature of 
factorization in QCD allows some of these assumptions 
to be justified by the A/ Em expansion. However, it 
should be noted that a priori there is no reason to prefer 
these factorization predictions to others that follow from 
the A/m& expansion (such as the prediction that certain 
strong phases are small). 

In Ref. [13 a x -fit was performed with 7 as a fit 
parameter, including decays to 77 and rf. The result 
7 = 61° ± 11° agrees well with global CKM fits. Here 
evidence for deviations from the standard model would 
show up as large contributions to the % 2 . The most recent 
analysis [l| has Br(K+n-), Br(K°ir°), and A CP (A' 7r ) 
contributing A\ 2 = 2.7, 5.9, and 2.9 respectively, giv- 
ing some hints f or p ossible deviations from the stan- 
dard model. Ref. [TtI extracted hadronic paramters from 
B — ► 7T7T decays, and used these results together with 
SU(3) and the neglect of exchange, penguin annihilation, 
and all electroweak penguin topologies except for the tree 
to make predictions for B — > Kir and B — > KK decays. 
They find large annihilation amplitudes, a large phase 
and magnitude for an amplitude ratio C/T which is in- 
terpreted as large P ut penguin amplitudes. The devi- 
ation of Br(Ar + 7r _ )/Br(A >0 7r°) from standard model ex- 
pectations was interpreted as evidence for new physics in 
electroweak penguins. 

There has been tremendous progress over the last few 
years in understanding charmless two-body, non-leptonic 
B decays in the heavy quark limit of QCD I25I l2(jL I27I 

E3, H3, 113, 13 HE E3, HE El IS M m- In 

this limit one can prove factorization theorems of the 
matrix elements describing the strong dynamics in the 
decay into simpler structures such as light cone distribu- 
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tion amplitudes of the mesons and matrix elements de- 
scribing a heavy to light transition |25| (for earlier work 
see Refs. 26]). It is very important that these results 
are obtained from a systematic expansion in powers of 
Aqcd/w&. The development of soft-collinear effective 
theory (SCET) [jj allowed these decays to be treated in 
the framework of effective theories, clarifying the separa- 
tion of scales in the problem, and allowing factorization 
to be generalized to all orders in a s . In Ref. [27j a proof 
of factorization was given for B — > DM~ type decays. 
Power corrections can also be investigated with SCET 
and in Ref. |2S| a factorization theorem was proven for 
the color-suppressed B° — > D°M° decays, and extended 
to isosinglet light mesons in Ref. j^H. Predictions from 
these results agree quite well with the available data, in 
particular the prediction of equal rates and strong phase 
shift for D and D* channels. 

Factorization for B — > M\ Mi decays involves three dis- 
tinct distance scales m 2 > E M A > A 2 . For B — > MiM 2 
decays a factorization theorem was proposed by Beneke, 
Buchalla, Neubert and Sachrajda [2j|, often referred to 
as the QCDF result in the literature. Another proposal is 
a factorization formula which depends on transverse mo- 
menta, which is referred to as PQCD 1^9. The factoriza- 
tion theorem derived using SCET |29|,|32j agrees with the 
structure of the QCDF proposal if perturbation theory 
is applied at the scales m 2 and mt,A. (QCDF treats the 
cc penguins perturbatively, while in our analysis they are 
left as a perturbative contribution plus an unfactorized 
large 0(v) term.) Due to the charm mass scale the iden- 
tification of a convergent ex pan sion for the cc penguins 
remains unclear [sTL IsM l34l l35j . For further discussion 
see [13, ) The SCET result improved the factoriza- 
tion formula by generalizing it to allow each of the scales 
m 2 , Em A, and Aq CD to be discussed independently. In 
particular, it was possible to show that a reduced set 
of universal parameters for these decays can already be 
defined after integrating out the scale m 2 [33, opening 
up the ability to make predictions for nonleptonic decays 
without requiring an expansion in a s (yEA). (If the m 2 
and m^A scales were separated in pQCD then this same 
result would be found for this first stage of factoriza- 
tion.) As a secondary step, additional predictions can 
be explored by doing a further expansion in a s at the 
intermediate scale. The expense of the second expansion 
comes in principle with the benefit of a further reduction 
in the number of hadronic parameters and additional uni- 
versality. In this paper we will explore the implications 
the first step of factorization has for B — > Ktt decays. 

There are several ways results from factorization can 
be used to analyze the data depending on i) whether 
perturbation theory is used at the intermediate -Em A 
scale as mentioned above, and ii) whether light-cone sum 
rules, models, or data is used to determine the hadronic 
parameters. In the QCDF 25] and PQCD [U analy- 
ses perturbation theory is used at the scale Em A- and 
light-cone sum rules |42l |43| or simple estimates were 
used for numerical values of most of the hadronic param- 



eters. Nonleptonic decay have also been studied with 
light-cone sum rules 0. With this input, all nonlep- 
tonic observables can be predicted and confronted with 
the experimental data. In both QCDF and PQCD a sub- 
set of power corrections are identified, parameterized in 
terms of new unknowns, and included in the numerical 
analysis. These power corrections are crucial to get rea- 
sonable agreement with the data. In these analyses it 
is sometimes difficult to distinguish between the model 
independent predictions from the heavy quark limit and 
the model dependent input from hadronic parameters. 
Ciucchini et al. have argued that so called charming 
penguins could be larger than expected and include un- 
knowns to parameterize these effects [3lJ. Fitting the 
hadronic parameters to non-leptonic data in some chan- 
nels and using the results to make predictions for other 
channels, as we advocate in this paper, has the advan- 
tage of avoiding model dependent input. Fits in QCDF 
have been performed in 10]. So far restrictions on the 
size of leading and subleading hadronic parameters nec- 
essary to guarantee convergence have not been explored. 
Other fits based purely on isospin symmetry have been 
explored in 0, [33] ■ 

In Ref. [32I l45l |4(| the factorization theorem was used 
in a different way, focusing on B — > rnr decays. Here 
perturbation theory was only used at the m 2 scale and 
fits to nonleptonic data were performed for the hadronic 
parameters in the LO factorization theorem. The prob- 
lematic contributions from charm-quark penguins were 
treated using only isospin symmetry. (This is also a 
good approach if power corrections spoil the expansion 
for this observable. Note that it avoids expanding the 
amplitude which has possible contamination from "chi- 
rally enhanced" power corrections [25|.) Here we con- 
tinue this program for B — > Kir and B — > KK decays 
(along with there comparison with B — > nn). For sim- 
plicity we refer to this as an "SCET" analysis, although it 
should be emphasized that other approaches to using the 
SCET-factorization theorem are possible. A key utility of 
factorization for nonleptonic decays is that the K/E and 
a s (mb) expansions are systematic and give us a method 
to estimate the theory uncertainty. Based on these un- 
certainties we investigate if the theory at leading order 
is able to explain the observed data. When deviations 
are found there are several possible explanations, all of 
which arc interesting: either the expansions inherent in 
the theoretical analysis are suspect, or there are statisti- 
cal fluctuations in the data, or we are seeing first hints 
of physics beyond the standard model. 

This paper is organized as follows: In section[H]we dis- 
cuss the theory input required to describe the decays of 
a B meson to two light pseudoscalar mesons. We briefly 
review the electroweak Hamiltonian at /1 = m& and then 
we discuss the counting of the number of parameters re- 
quired to describe these decays using SU(2), SU(3), and 
SCET analyses. We finish this section by giving a gen- 
eral parameterization of the decay amplitudes in SU(2). 
(In the appendix we give the relations between our pa- 
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rameters and the graphical amplitudes |23j, |47J, |48(.) In 
section lllll we give the expressions of the decay ampli- 
tudes in SCET. We begin by giving the general expres- 
sions at leading order in the power expansion, but correct 
to all orders in a s and comment about new information 
that arises from combining these SCET relations with the 
SU(3) flavor symmetry. We then use the Wilson coeffi- 
cients at leading order in a s (m b ) and give expressions for 
the decay amplitudes at that order. We finish this sec- 
tion with a discussion of our estimate of the uncertainties 
which arise from unknown 0(a s (m b )) and 0(Aqcd/E) 
corrections. A detailed discussion of the implications of 
the SCET results is given in section [fVl We emphasize 
that within factorization the ratios of color suppressed 
and color allowed amplitudes (C/T and EW C / EW T ) 
can naturally be of order unity at LO in the power count- 
ing, contrary to conventional wisdom |3^ |. We also per- 
form an error analysis for the Lipkin and CP-sum rules in 
B — > Kir decays, and discuss predictions for the relative 
signs of the CP asymmetries. We then review the infor- 
mation one can obtain from only the decays B — » 7T7t, 
before we discuss in detail the implications of the SCET 
analysis for the decays B — > K ir and B — > KK. 



II. THEORY INPUT 

A. The electroweak Hamiltonian 

The electroweak Hamiltonian describing Ab = 1 tran- 
sitions b — ► / is given by 
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where the CKM factor is X ( ,P = V pb V* f . The standard 
basis of operators are (with 0\ <-> 0\ relative to p9|) 

0\ = (pb)v-A(fp)v-A, P 2 = (ppb a )v-A(f a Pp)v-A, 

O3.4 = {(pb) V -A(qq)v-A , (7/3 b a ) V -A{q a qp)v-A}, 

05,6 = {Q b )v-A{qq)v+A , (7 pb a ) V - A (q a qp)v+A} , 

3e — — 

07,8 = -Yi(fb)v-A(qq)v+A,(fpb a ) V -A(q a qp)v+A}, 

3e — — 

°9,io = -^{{fb)v-A{qq)v-A , (//c,)v-a(m)i'-a}, 

°77,s 9 = J^{eF^gG^T a }(l+ 75 )b. (6) 

Here the sum over q — u,d,s,c,b is implicit, a, (3 are 
color indices and e q are electric charges. The AS = 
and AS = 1 effective Hamiltonian is obtained by setting 
f = d and / = s in Eqs. (TiHil) . respectively. The Wilson 
coefficients are known to NLL order ^^j. At LL order 
taking a s (m z ) = 0.118, m t = 174.3, and m b = 4.8 GeV 
gives C7 7 (TOb) = —0.316, C% g {mb) = —0.149 and 

Ci_i (to 6 ) = {1.107 , -.249 , .011 , -.026 , .008 , -.031 , 
4.9 x 10~ 4 , 4.6 x 10~ 4 , -9.8 x 10~ 3 , 1.9 x 10~ 3 } . (7) 



Below the scale /i ~ one can integrate out the bb 
pairs in the operators O3-10 The remaining operators 
have only one 6-quark field, and sums over light quarks 
q = u,d, s, c. This gives rise to a threshold correction to 
the Wilson coefficients, 



C~(m b )=C+(m b ) 



1 



a s {m b ) 

47T 



Sri 
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+ ~r 6r c 

47T 



(8) 



where C + and C are the Wilson coefficients with and 
without dynamical b quarks, and SrJ and 5r^ are given 
in Eqs. (VII.31) and (VII.32) of 49]. This changes the 
numerical values of the Wilson coefficients by less than 
2%. Integrating out dynamical b quarks allows for addi- 
tional simplifications for the electroweak penguin opera- 
tors, since now for the flavor structure we have 
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^e q (fb)(qq) 



- [fb)(2uu -dd- ss + 2cc) 



(9) 



\{Jb)(cc)- 1 - £ (fb)(qq). 
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q—u,d 7 s,c 



The operators Og and O10 have the regular (V — A) x 
(V — A) Dirac structure, and can therefore be written as 
linear combinations of the operators O1-4, 



Oi0 = \oi + ~Ol- l -O A 



(10) 



This is not possible for the operators O7 and Os, which 
have (V—A) x (V+A) Dirac structure. Thus, integrating 
out the dynamical b quarks removes two operators from 
the basis. To completely integrate out the dynamics at 
the scale m b we must match onto operators in SCET, as 
discussed in section |Hl] below. 



B. Counting of Parameters 

Without any theoretical input, there are 4 real 
hadronic parameters for each decay mode (one complex 
amplitude for each CKM structure) minus one overall 
strong phase. In addition, there are the weak CP vio- 
lating phases that we want to determine. For B — ► ttt: 
decays there are a total of 11 hadronic parameters, while 
in B — > Kit decays there are 15 hadronic parameters. 

Using isospin, the number of parameters is reduced. 
Isospin gives one amplitude relation for both the tttt 
and the Kit system, thus eliminating 4 hadronic pa- 
rameters in each system (two complex amplitudes for 
each CKM structure). This leaves 7 hadronic parame- 
ters for B — > 7T7T and 11 for B — > Kir. An alternative 
way to count the number of parameters is to construct 
the reduced matrix elements in SU(2). The electroweak 
Hamiltonian mediating the decays B — > tttt has up to 
three light up or down quarks. Thus, the operator is 
either AI — 1/2 or AI — 3/2. The two pions are ei- 
ther in an / = or I = 2 state (the 1 = 1 state is 
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no 
expn. 


SU(2) 


SU(3) 


SCET 

+SU(2) 


SCET 
+SU(3) 


B — > 7T7T 


11 


7/5 


15/13 


4 


4 


— > Ktt 


15 


11 


+5(6) 


B -> KK 


11 


11 


+4/0 


+3(4) 


+0 



TABLE II: Number of real hadronic parameters from differ- 
ent expansions in QCD. The first column shows the number of 
theory inputs with no approximations, while the next columns 
show the number of parameters using only SU(2), using only 
SU(3), using SU(2) and SCET, and using SU(3) with SCET. 
For the cases with two numbers, #/#, the second follows 
from the first after neglecting the small penguin coefficients, 
ie setting C 7 ,s = 0. In SU(2) + SCET B — > Ktt has 6 pa- 
rameters, but 1 appears already in B — > hit, hence the +5(6). 
The notation is analogous for the +3(4) for B —* KK. 



ruled out by Bose symmetry). This leaves 2 reduced 
matrix elements for each CKM structure, (0||l/2||l/2) 
and (2||3/2||l/2). For B -> Ktt decays the electroweak 
Hamiltonian has either AI = or AI = 1. The Ktt 
system is either in an / = 1/2 or I — 3/2 state thus 
there are three reduced matrix elements per CKM struc- 
ture, (3/2||l||l/2), (l/2||l||l/2) and (l/2||0||l/2). Fi- 
nally, KK is either an / = or / = 1, and there are 
again three reduced matrix elements per CKM structure, 
<0||l/2||l/2), (l||l/2||l/2),and (l||3/2||l/2). 

The SU(3) flavor symmetry relates not only the decays 
B — > 7T7T and B — > Ktt, B — > KK, but also B — > Trrjg, 
B — >• rj 8 K and B s decays to two light mesons. The de- 
composition of the amplitudes in terms of SU(3) reduced 
matrix elements can be obtained from |50L l5ll |52j . The 
Hamiltonian can transform either as a 3 , 3 , 6 or 15. 
Thus, there are 7 reduced matrix elements per CKM 
structure, (1||3 S ||3), (1||3 Q ||3)^(8||3 S ||3), j8||3~ a ||3), 
(8||6 S ||3), (8||15 ||3) and (27||15 s ||3). The 3° and 3 s 
come in a single linear combination so this leaves 20 
hadronic parameters to describe all these decays minus 1 
overall phase (plus additional parameters for singlets and 
mixing to properly describe r\ and 77'). Of these hadronic 
parameters, only 15 are required to describe B — ► -ktt 
and B — » Kit decays (16 minus an overall phase). If we 
add B — > KK decays then 4 more paramaters are needed 
(which are solely due to electroweak penguins). This is 
discussed further in section ITTDl 

The number of parameters that occur at leading order 
in different expansions of QCD are summarized in Ta- 
ble HH including the SCET expansion. Here by SCET 
we mean after factorization at rrib but without using any 
information about the factorization at y/EA. The SCET 
results are discussed further in section ITTT1 but we sum- 
marize them here. The parameters with isospin+SCET 
are 

tttt: {C^K^/Uf^™}, (11) 
Ktt : {C^ + C? \ 0K<f", C BR + C?*, 0«C? R ,Pk«} , 



Here PmiM 2 are complex penguin amplitudes and the re- 
maining parameters are real. 1 In B — > 7T7t the moment 
parameter p„ is not linearly independent from the pa- 
rameters ( B7r and C, B 77 , and only the product P-kCj^ was 
counted as a parameter. In any case it is fairly well known 
from fits to 7*7 -► tt° [H3 3/3 w = {x- 1 )* ~ 3.2 ± 0.2. In 
isospin + SCET B — > Ktt has 6 parameters, but the first 
one listed in (|llfl appears already in B — > 7T7t, hence the 
+5 in Table ITT1 If the ratio (3k /P* was known from else- 
where then one more parameter can be removed for Ktt 
(leaving +4). For B — > KK we have 4 SCET parame- 
ters. One of these appears already in B — ► Ktt, hence 
the +3, and if Pk/Pk is known from other processes it 
would become +2. 

Taking SCET + SU(3) we have the additional relations 

fB-K _ aBK _ /-BK sBi: _ /-BK _ rBK a _ a _ a _ 

and = = A^ C K which reduces the number of 

parameters considerably. 

Note that there are good indications that the param- 
eters ( BM and Cj M are positive numbers in the SCET 
factorization theorem. (Pk, P-k, Pk are a l so positive.) 
This follows from: i) the fact that are related 

to form factors for heavy-to-light transitions which with 
a suitable phase convention one expects are positive for 
all q 2 , ii) that £f M is positive (from the relatively safe as- 
sumption that radiative corrections at the scale V EA do 
not change the sign of ^iM 2 an( ^ q ^ p^\ B > 

and finally iii) that the fit to B — > 7T7r data gives 
(•Stt^Stt > so that SU(3) implies ( BK ,(j K > 0. We 
will see that this allows some interesting predictions to 
be made even without knowing the exact values of the 
parameters. 

In using the expansions in © it is important to keep 
in mind the hierarchy of CKM elements, and the rough 
hierarchy of the Wilson coefficients 



Ci > C2 3> C3_6 ^> Cg 10 <L C*7 : g 



(12) 



Some authors attempt to exploit the numerical values of 
the Wilson coefficients in the electroweak Hamiltonian 
to further reduce the number of parameters. A common 
example is the neglect of the coefficients CV^ relative to 
Cg^o- In Eq. (|TU1) the electroweak penguin operators O9 
and O\o were written as linear combinations of O1-4. 
This implies that if one neglects the electroweak penguin 
operators Q7 and Qs, then no new operators are required 
to describe the EW penguin effects. In some cases this 
leads to additional simplifications. One can show that for 
B — > 7T7T decays the A J = 3/2 amplitudes multiplying the 
CKM structures X u and A c are identical [23,|24(. Thus, 



KK-. {t BK +cr^Ktr,p K K} 



BK 



-BK 



1 The penguin amplitudes are kept to all orders in A/ mi, since so 
far there is no proof that the charm mass m c does not spoil fac- 
torization, with large a s (2m c )v contributions competing with 
a s ( mi,) ha rd-charm loop corrections |32| . This is controver- 
sial [2,4, 35]. Our analysis treats these contributions in the most 
conservative possible manner. 
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SU(2) gives one additional relation between complex am- 
plitudes in the tttt system, reducing the hadronic param- 
eters to 5. For B — > Ktt decays the operators giving 
rise to the A 3 / 2 reduced matrix elements are identical for 
the X u and A c CKM structures only if SU(3) flavor sym- 
metry is used [Hij . Thus, for these decays two hadronic 
parameters can be eliminated after using SU(3), leaving 
13. Considering B — > KK adds two additional parame- 
ters. Note that dropping C^s makes it impossible to fit 
for new physics in these coefficients. In our SCET analy- 
sis all contributions C7-10 are included without needing 
additional hadronic parameters. 

Finally, some analyses use additional "dynamical as- 
sumptions" and drop certain combinations of reduced 
matrix elements in SU(3). For example, the number 
of parameters is often reduced by neglecting parameters 
corresponding to the so called annihilation and exchange 
contributions. 



C. General parameterization of the amplitudes 
using SU(2) 

Using the SU(2) flavor symmetry, the most general am- 
plitude parameterization for the decay B — > tttt is 

A(B° - iT+n-) = -AWT™ - AMP„ (13) 
-» ttV) -Af>C„ -X^(EW^ - P ww ) 
V2A(B- -» ir-ir°) = -X ( u d) (T^ + CW) - \ { c d) EW^ 

where we have used the unitarity of the CKM matrix 
X { t f) = -A ( / } - \ { J ] . The amplitude parameter EW T ^ 
receives contributions only through the electroweak pen- 
guin operators 0-?-xo- For B — > Ktt decays we write 

A(B-^n-K Q ) = X^A K7T + X^P Kn (14) 
V2A(B- -> tt°K-) = -\i s) (C K ir+T K7r + A K7r ) 
-X^(P K7r + EW Kw ) 
A(B ^7t+K-) = -X^T Kv 

-X^(P Kv + EW K7r ) 
V2A(B° - 7T°K°) = -X^Ck* 

+Xi s \P Kn -EW Kw + EW K% ) 

Finally for B — > KK decays there is no SU(2) relation 
between the amplitudes and we define 

A(B-^K-K°) = X^A KK + X^P KK 
A(B° — > K a K a ) = X^B KK 

+X { c d) (Pkk + PAkk + EWkk) 
A(B° -> K-K+) = X^E KK - X^PA KK . (15) 

As mentioned before, after eliminating X\ there are 
four complex hadronic parameters for B — > tttt and six 
for B — > Ktt. The additional relation one obtains in the 



limit C7 8 — > is 



EWt 



3 C 9 + C 



in 



yP'K'K ~t~ C^^) ; 



(16) 



2 C x + C 2 

where we have neglected terms quadratic in Cg or C\q. 

The EW amplitudes are purely from electroweak pen- 
guins, however there are also electroweak penguin con- 
tributions in the other amplitudes as discussed further in 
section lTlI Dl Also, the hadronic parameters in Eqs. i|13|) - 
(|15JI are a minimal basis of isospin amplitudes, not graph- 
ical amplitude parameters. In the appendix we show how 
these amplitude parameters are related to the graphical 
amplitudes discussed in *23l I47L lisj . 



D. Additional relations in the SU(3) limit 

In the limit of exact SU(3) flavor symmetry the param- 
eters in the tttt system and the Ktt system satisfy the two 
simple relations [47 * 



EWl v = EW T K « 



(17) 

Thus, the hadronic parameters in the combined Ktt, tttt 
system can be described by 8 complex parameters (15 
real parameters after removing an overall phase), if no 
additional assumptions are made. A choice for these pa- 
rameters is 



T, 
EW' 



TTTT 1 

T 

K-rr > 



EW 



c 

Km 



P A 

1 7T7r ) 

AC, 



Ktt 1 

AP 



where we have defined 

AC = Ck-tt ~ Crjr 

AP = P^ - P K7r 



(18) 



(19) 



This can also be seen by relating these amplitude param- 
eters directly to reduced matrix elements in SU(3), which 
can be done with the help of the results in Ref. ■ As 
before, if the small Wilson coefficients C7 and C$ are 
neglected, we can again use the relation in Eq. (|TrJ)l to 
eliminate one of the 8 complex hadronic parameters. 
Four additional relations exist if the amplitudes for 



B — > KK are included 



A K k 


— Aktt 




Pkk 


= Pkk 




Ekk 


= Tktx~ 


-T — 

7T7T 


PAkk 


— p - 

J 7T7T 


Pk.- 



(20) 



-AC 



AP 



EW^ 



In the limit of vanishing Wilson coefficients Cj and C% 
there are two additional relation 23] 
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C 9 + C 



10 



Ci 



EWkk 



f C 2 
3C^ 
2 Cx 



Cx — C-2 

{Aktt - T, 



C 2 



(Ak-k — Ptttt + Cktt + Pkk) 
Cktt — ZCtttt — Pkk) 
{B kk - A KK + Ekk) (21 
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E. Sum-Rules in B — ► Kn 

In this section we review the derivation of two sum- 
rules for B — ► ifor, the Lipkin sum-rule [H^, |HE E3| and 
CP sum-rule j^. Higher order terms are kept and will 
be used later on in assessing the size of hadronic correc- 
tions to these sum rules using factorization. To begin we 
rewrite the SU(2) parameterization of the amplitudes as 



A(B~ 
A(B° 
V2A(B~ 



= \ { ° ] Pk. 

-> TT + K~) 



(22) 



1 - - eA e-^e*^ 



= -AWPj fw [l + -(eg B e < *™ 



-ere 



pT—vy) 



V2A(B° - 7T°K ) 



1+ ffie* 



-ee 



20 — 27 > 



1 ^etu ^ 



20 e 



-ece 



"PC— '71 



where 



2 ere 



2 £ce 



i<pA 



A« 


( 


-Tktt) 


A« 




Pk, ' 


An 


( 


-Cktt) 






Pk* ' 


\W 
a« 


( 


-A Kn ) 






Pktt 


\W 

An 


( 


-Tk-k — 


A« 





(23) 



Pk-k 



and 



1 _eto i</ 

2 £t 6 

- e5T e l<# 

2 c 

1 

2 ^eu> ^ 



A'tt 



Ptf,r 



(24) 



At: 



i0 e 



Pk,t ' 

A 71 



At, 



Pktt 



These parameters satisfy 



ew %4 



~ e C 



(25) 



The non-electroweak e-parameters are suppressed by the 
small ratio of CKM factors /\ { c ] \ ^ 0.024 but are 
then enhanced by a factor of ~ 4-15 by the ratio of 
hadronic amplitudes. The electroweak e-parameters are 
simply suppressed by their small Wilson coefficients and 
end up being similar in size to the non-electroweak e's. 



Next we define deviation parameters for the branching 
ratios 



Ri 
Ra 



2Br(5~ 


-* TT°K-) 


Br(B- - 


-> TT-K°) 


Br(B° - 


y tt-K+)t b - 


Bt(B~ - 


-> tt-K°)t b o 


2Br(S° - 


-> ir°K )r B - 


Br(B- - 


-» tt-K°)t B o 



(26) 



-1, 
-1, 



(27) 



and also rescaled asymmetries 

Ai = {l + R 1 )A CP (n°K-), 
A 2 = {l + R 2 )A CP (jr-K+), 
A 3 = (l + R 3 )A CP (ir°K ), 
A 4 = A CP (n-K°). 



The division by Bi^tt - ^ ) in the A$ asymmetries is not 
necessary but we find it convenient for setting the nor- 
malization. Expanding in ca we find that to second order 
in the e parameters the Ri are 

Ri = [e™ cos <j>™ — e cos cos 7 + ca cos 4>a cos 7] 

' ' " ^t 2 ~^a)-\^ c t 'cos7cos(0-0^) 



R-2 



Ra 



+ (ef^ cos 0^*" — e cos cos 7) ca cos a cos 7 
+ e 2 A cos 2 0a cos 2 7 , 

[eg" cos 0^ — ct cos 0t cos 7+ £a cos 0a cos 7] 

" 1 (4 + e - 2 _ £ 2 ) _ 1 £t e ej CQS 7 cos( ^ _ ^ 

+ (eg" cos 0^ — er cos 0t cos 7) ea cos a cos 7 
+ e A cos 2 4>a cos 2 7 , 

[ — e ew cos ew + ec cos c cos 7 + €a cos ^a cos 7] 
(e ew 2 +e e J 2 -e A )-^e c cos 7 cos(0 etu - C 
(e e ™ cos (j) ew — ec cos 0c cos 7) ca cos a cos 7 



e A cos 2 4>a cos 2 7 



(28) 



and the A, are 



Ai = sin7[ — esii 



sm(0-0™)] 



[l + €A COS 0A COS 7], 



sin 7 [ — ex sin 0^ — — ex e & c sin(0x - 
x [l + ea cos (f>A cos 7] j 

. r . , 1 



sin 7 |ec sm 0c - 2 ec e 6 '"' si 1 ) 



X [l + €A COS (j)A COS 7] , 

sin 7 [ — ca sin 0a] [l + £A cos 0a cos 7] . (29) 
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Note that these rescaled CP- Asymmetries are indepen- 
dent of the electroweak penguin e's at 0(e). This is not 
true for the original asymmetries Aqp- 

Sum rules are derived by taking combinations of the 
Ri and A< which cancel the 0(e) terms. The Lipkin sum 
rule is the statement that 

Ri -R2 + R3 =0(e 2 ) 

--( e 2 e 2 + e 2 f 2 ^ e ew2 f ew2 , ew2\ 

— - e T + £ C ~ e A + £ T - e C + e ) 

+ e^ cos 2 (0a) cos 2 (7)- - cos(7)[e^ecos(^-^) 

- e T e e c w cos(0 T -</>§") + e c e ew cos(0 c -^)] , (30) 

where we used the real part of Eq. (|25ll . The CP-sum 
rule is the statement that using the imaginary part of 
Eq. I)23H the 0(e) terms cancel in the sum 

A 1 -A 2 +A 3 -A 4 = 0(e 2 ) 

= -1 sin( 7 ) [e^esin(0 - - e T e e J> sin(0 T - ^c) 



e c e ew sin(4>c - 4>ew)] 



(31) 



The accuracy of these sum rules can be improved if 
we can determine these 0(e 2 ) terms using factorization. 
This is done in section llVD II 



III. AMPLITUDE PARAMETERS IN SCET 



A. General LO expressions 

The factorization of a generic amplitude describing the 
decay of a B meson to two light mesons, B — > M1M2, 
has been analyzed using SCET [32|. Here Mi and M 2 
are light (non-isosinglet) pseudoscalar or vector mesons. 
The SCET analysis involves two stages of factorization, 
first between the scales {nib or Em} 2 ^> -EmAqcd, and 
second between £?mAqcd 3> Aq Cd . Here we only con- 
sider the first stage of factorization where we integrate 
out the scales {nib, Em}, and keep the most general pa- 
rameterization for physics at lower scales. It was shown 
in Ref. [32( that a significant universality is already ob- 
tained after this first stage, in particular there is only one 
jet function which also appears in semileptonic decays to 
pseudoscalars and longitudinal vectors. This leads to the 
universality of the function we call C, BM (z). We note 
that this also proves that the second stage of matching is 
identical to that for the form factor, so the SCET results 
for form factors in Refs. 59!] can immediately be applied 
to nonleptonic decays if desired. A summary of the anal- 
ysis of SCET operators and matrix elements is given in 
Appendix [X] 

After factorization at the scale mb the general LO am- 



plitude for any B — > M1M2 process can be written 
GFm 2 B 



A 



V2 



{fMiJdu dz T u (u, z)Cj Ma {z)4> M > (u) 



+f Ml C M2 duT- t 
Jo 



' K {u)^{u)\ + {\^2) 



(32) 



where ( BM and Cj M are non-perturbative parameters 
describing B — ► M transition matrix elements, and 
A^. lM2 parameterizes complex amplitudes from charm 
quark contractions for which factorization has not been 
proven. Power counting implies C, BM ~ Cj M ~ 
(A/Q) 3 / 2 . Tu{u,z) and T\q(u) are perturbatively cal- 
culable in an expansion in a s (mb) and depend upon the 
process of interest. 

It is useful to define dimensionless hatted amplitudes 



i = 4(GeV" 1 ) 



No 



Gpm 2 B 

~7T 



(33) 



Using Eq. (|32J) we find that the amplitude parameters in 
the B — > 7T7T system are 



^7T 



ew 
c lt 



-fn [(Clu + C4u 
~U [(C2u 



-2t 



Lew 



C4t — C4u)ttC 

+ ((b2u-b c 2 y + b™ + b 4t -b 4u )C 
-A™ + U[(c<% + c it U B * 

f \ lrP w -X- r- cw — r cw \ f Brr 
Jt~ L\ It 2t L 3t /rt 

+((bl7 + bl7-btT)t? v U • 



Btv 



?*)*] 



(34) 



For the B — » Ktt system we find 



TK n — —fx [(Clu 

Cktt 



It - c 4t + Ciu) K (, 

+((b lu -btT-h t + b iu ){f«) R ] 



-f* 



ew\ /-BK 



(C2u - c c 2 r + c c 3 r).c 



Pk* 

Akt\ 



+{(b2u - &§T + btJ)cf% 

Ik [(at - C4u)/K s " + ({bu - MCf") 
: -A£< + Ik [(cu)rC B " + (bittf Ik] 
■ Ik [(dt - c 4u ) R C B7: + ((bit - MC?")i?j 



K 



£ f / CW\ f-B'7T 1 /LCW /-BlT\ 

Ik [(c u } k C + (bu Cj )k\ 



1 f /„cw _ ew\ /-BK 
~r Jit \<- 2t t- 3t /w^ 

+ ((bTt-bTt)Cj 

EW c K ^ = i K [(ciY) R c B7r + (bl7(:f lK 



.7 /w 



(35) 
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Note that the dominant non-factorizable charm elec- 
troweak penguin contribution, is absorbed into A cc . Fi- 
nally, for the B — > KK system we find 



Akk =B 

Pkk 
Erk 



KK 



(36) 



K 



(c4t - cau)kC BK + ({bit - bi u ) Cj K )k 



BK\ 



_ A KK 

PA KK 



+ fK (c 4t ) K C BR + (b it Cj R )K 



EW 



KK 



0. 



In Eqs. (|34l3t)|) we have defined 



(Ci)M = I du Ci(u)(f)M(u) 

Jq 



(37) 



du / dz b i (u,z)cf>M 1 (u)(j 2 {z). 
o Jo 



We have also decomposed the Wilson coefficients of 
SCET operators defined in [33 as 



(38) 



and in some equations we have split the contributions 
from strong (0 3 -e) and electroweak penguin operators 
(O7-10) to the cu and bu 



b lt = bi+b™. 



(39) 



Note that to all orders in perturbation theory one has [48 



c lt — c 2t 

b p lt = b p 2t - 



= c 3u 



= 4t 

b P 3t = 



0. 



(40) 



All the hadronic information is contained in the SCET 
matrix elements C S,r , Cj 7r ) C BK ', Cf K , and 
the decay constants / ff and fx, and the light cone distri- 
bution functions of the light mesons cf> n (x), 4>k{%), and 
4>k(x). 

The Wilson coefficients and 6^ are insensitive to the 
long distance dynamics and can therefore be calculated 
using QCD perturbation theory in terms of the coeffi- 
cients of the electroweak Hamiltonian, Cj. Any physics 
beyond the standard model which does not induce new 
operators in H\y at /i = mw will only modify the values 
of these Wilson coefficients, while keeping the expressions 
for the amplitude parameters in Eqs. (|34H36I) the same. 

We caution that although the amplitudes Ak-k and 
Akk do get penguin contributions at this order, they will 
have subleading power contributions from operators with 
large Wilson coefficients that can compete. Therefore 
their leading order expressions presented here should not 
be used for numerical predictions. As mentioned earlier, 
the penguin amplitudes are kept to all orders in A/ E. In 
all other amplitudes the power corrections are expected 



to be genuinely down by A/E when the hadronic pa- 
rameters are of generic size. In the observables explored 
numerically below it will be valid within our uncertain- 
ties to drop the small A^tt and Akk amplitudes and so 
this point will not hinder us. 



B. SU(3) limit in SCET 

In the SU(3) limit the hadronic parameters for pions 
and kaons are equal. This implies that 



r Bn _ f-BK 
SJ — W j 



qBtt _ ^BK 

{Ci)K = (Ci) w , (41) 
Furthermore 



A Kir _ A1TTT _ A 



KK 



(42) 



To see this note that in SCET the light quark in the oper- 
ator with two charm quarks is collinear and can therefore 
not be connected initial B meson without further power 
suppression. Without the use of SCET this so called 
"penguin annihilation" contribution would spoil the re- 
lation in Eq. {35)1. 

Using this we find two additional relations which are 
not true in a general SU(3) analysis but are true in the 
combined SCET + SU(3) limit 

AC = Or* - CW = , (43) 
AP-EW Kn = P ff7r - P K „ - EW Kn = , 

where the zeroes on the RHS are 0(m s /A) + 0(A/E). 
Using the SU(3) relation in Eq. I|20|l we see that these 
amplitudes are equal to "exchange" or "penguin annihi- 
lation" amplitudes that are power suppressed in SCET. 



C. Results at LO in a s (mb) 

While the Ci are known at order a s , the bi are cur- 
rently only known at tree level. For consistency, we thus 
keep only the tree level contributions to the Ci as well. In 
this case they are independent of the light cone fraction 
u and thus Ci(u) = c,, and there occurs a single nontriv- 
ial moment of the light-cone distribution function from 
the bi terms. Since the parameter A^[ lM2 cx a s (2m c ) it 
would be inconsistent to drop the a s corrections in the 
penguin amplitudes. However, as long as we have the free 
complex parameter A^f c lM2 these corrections are simply 
absorbed when we work with the full penguin amplitudes 
-P/W1M2 using only isospin symmetry. This is also true of 
chirally enhanced power corrections in PmiM 2 - 

Using LL values for the Wilson coefficients we find for 
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the non-electroweak amplitudes at fi — rut = 4.8 GeV 
Ci 



C2v 



Ci + ^ = 1.025 
C 2 + ^ = 0.121 



(44) 



4 = -(C 4 + ^) +0(da s ) = 0.022 + G(a s ). 

Here 0(Cia s ) indicate unsuppressed a s corrections that 
were computed in Ref. [2^| and verified in |2|j . The con- 
tributions from the operators O7-10 give 



-2( 



-3t 



ow 
' 4/ 



' s (r + ° 9 

"2r 9+ AT 

2{ Cl0 + K 



= 0.0021 

= 0.0138 
= -0.0010 

-0.00068 . (45) 



The coefficients bi(u, z) are independent of the variable 
z at leading order and we write bi(u, z) = bi(u). For the 
non-electroweak amplitudes we have 



hu{u) 



b 2 u(u) 



= 1.025 



u) N c 
0.249 



(46) 



3u 



c 2 + (i + 1)£l 



0.121 



1.107 

3m 



b\ t {u) = -C4-(l + ^)| + 0(Cio.) 

= 0.022- -£H + 0(a s ) 

where C(Cia s ) denotes unknown unsuppressed a s cor- 
rections and 22 = 1 — it. For the clcctroweak terms 



blJ(u) 



b™{u) 



C 



10 



i\ c 9 
1 + - — 

22/ 7V r 



0.0021 4 



0.0147 

322 



10 



22/ N n 



0.0138 



0.0029 



6|T(«) 



-0.00010 

5 [Cio + I I 
-0.00068 



3m 

l_\Cs 
u) N r 
0.00069 
3u 

uJ N c 
0.0049 



3 11 



(47) 



Note that only &|™ involves u and that this coefficient 
only appears convoluted with pions in Eqs. ({34136(1 . For 
pions one can take 1/m — > 1/22 using charge conjugation 
and isospin. Since the 6's then only involve factors of 1/22 
it is useful to define the nonperturbative parameters 



Jo 



3u 



(48) 



Using these values for the Wilson coefficients we ob- 
tain the amplitude parameters in terms of the non- 
perturbative parameters in the tttt system at fi = mj, 

T ff7r = -0.131 (C S " + Cf") +0.031 0„tf* 
= -0.017(C S " + C.D -0.144 faCj* 



EW„ 
ft 



0.0022 (C Brr + (j n ) + 0.0015 
-AH + 0-0030 (C Brr + Cj w ) 
-0.0002 ArCj' + A^, (49) 



where A^ is the additional perturbative correction from 
a4t and &4t at 0(a s (rrib)) which can involve larger Wilson 
coefficients like C\$. (We could also include large power 
corrections in A^ assuming that such a subset could be 
uniquely identified in a proper limit of QCD.) We do not 
need knowledge of A B V for our analysis since there are 
two unknowns in each of P W7r and A™ and we will simply 
fit for Ptttt. 

In the Kir system we find 



Crtt 



-0.160(C BTr + Cj w ) + 0.040 (3 K (f * (50) 
-0.003(C BTr + £f*) + 0.003 PkCj* 
-0.014(C ba ' + Cj K ) -0.146 /3^C. 



EW Kn = 0.0019 (C BK + Cj K ) ~ 0.0005/3,rCj K 
+0.0003 (( B7! + ( B7r ) + 0.0023 PkCj* 

EW C Kl , = 0.0003(C s,r + Cf^) + 0.0023 /foCj* 
Arw = -A^ + 0.0034(C B7r +Cf r ) 



-0.0026 /3 K Cj + A 



'if* 



i K7r = 0.0034(C B,r + Cj 37r )-Q-0026feC 



where A^ and A^ are analogous corrections to A^. 
For Pr-,,- the perturbative correction competes with A^ . 
For Aktt power corrections could be in excess of the lead- 
ing order value, so that any numerical value for this am- 
plitude is completely unreliable at the order we are work- 
ing. 

Finally for the B — > KK amplitudes that have a con- 
tribution from the LO factorization theorem we have 



A KK = 0.0034(C BA ' + Cf") - 0.0026/^Cj 



BK\ 



-BK 



0(a 3 ) 



KK 



-A 



KK 



0.0034(C 



BK 



+ Cf K ) 



0.0026/3rtCj K + A 



KK 



(51) 
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and Bkk = Akk- Here the value of Akk is not reliable, 
since it will have large 0{C\ A/m*) power corrections that 
are likely to dominate. 



D. SCET Relations for EW penguin amplitudes 

Using the SCET results in the previous two sections 
it is simple to derive relations that give the electroweak 
penguin contributions in terms of tree amplitudes 



where dropping C3.4 relative to Ci 2 



1 M 1 M 2 



T, 



M 1 M 2 \ Cvio=0 , 



C M l M 2 - C MiM 2 lc 7 _ 10 =0 ' 



(52) 



Such relations are useful if one wishes to explore new 
physics scenarios that modify the electroweak penguin 
parameters C7_io in H\y- To separate out all electroweak 
penguin contributions in B — > mr and using SCET to- 

- c 

gether with isospin we define EW nv by 

c 

7T7T 5 



f^=f°+E~W 



P = P° 

± 7T7T tt' 



7T7T 1 

- C 

, EW . 

7T7T 1 7T7T 



C' 



At LO in SCET we find 

c 



e 2 CkTT 1 



EW^ = ei T°* 
Ewl v = e 3 f^ + e A C ^, 
where dropping relative to C\ t 2 one finds 
C10C1 — CgC2 



c-'i 



ci-ci 



-2.9x10" 



(53) 



(54) 



(55) 



C9C1— C10C2 m-3 

e 2 = — — — — = -8.9 X 10 



e3 



e 4 



Ci (3Cio - 3C 7 - 2C 8 + 3C 9 ) - 3C 2 (C 10 + C 8 + C 9 ) 



2(C 2 -Cf) 

-1.5 x 10~ 2 , 

- C 2 (3Ciq - 3C 7 - 2C 8 + 3C 9 ) + 3d (C 10 + C 8 + C g ) 
2(C 2 -Cf) 
= -1.3 x 10~ 2 . 

The numbers quoted here are for the standard model LL 
coefficients. 

For B —> Kn we separate out the electroweak penguin 
contributions by writing 

f Kn = f K „ + ^EWZ„ , A Kv = P K „ - \eW° k „ , 

C Kn = C% n + EW T K ^EW C Kn , 

P Kn = P^-±EW% n , (56) 

and find that SCET+isospin gives 

EWL = j; (e 5 + e 6 6°J) + e 7 f K „ , 



es 



e6 = 



3Ci(CiCg — C2C10) 

2C 2 (C 2 -C 2 ) 
3(CiqC2 — C9C1) 



= -6.0 x 10" 



2(C 2 



-ci) 



= -1.3 x 10" 



e 7 = || = 5.9xl0- 2 . 



(57) 



For the amplitude EW^-^ no such relation exists, if the 
inverse moments (3^ and (3k are taken as unknowns. One 
can still use the SU(3) relation in Eq.Elto equate EW T 
in the Kir and 7T7t system. 



E. Estimate of Uncertainties 

These expressions of the amplitude parameters are cor- 
rect at leading order in AQCD/'-E'-n-, and as we explained 
above, the complete set of Wilson coefficients is currently 
only available at tree level. Thus, any amplitude cal- 
culated from these SCET predictions has corrections at 
order a s (rrib) and Aqcd/^tt- Using simple arguments 
based on dimensional analysis, we therefore expect cor- 
rections to any of these relat ions at the 20% level. We are 
working to all orders in a s (v / AB), and so we avoid adding 
additional uncertainty from expanding at this scale. 

Note that we have allowed for a general amplitude 
Ptttt , which contributes to the reduced isospin matrix ele- 
ment (l/2||0||l/2) in the Kit system, and to the reduced 
isospin matrix element (0||l/2||l/2) in the 7T7t system. 
All power correction contributing to the same reduced 
matrix element will be absorbed into the value of the ob- 
servable P W7r . In the following will thus fit directly for the 
parameters and Pktt, which reduces the theoretical 
uncertainties significantly. This implies that the theoreti- 
cal uncertainties on the amplitude parameters P^, Pktt, 
and Pkk are ~ 3% from isospin rather than ~ 20%. All 
other appreciable LO amplitude parameters are consid- 
ered to have uncertainties at the 20% level. 

Using this information, we can now estimate the size 
of corrections to the individual observables. For the de- 
cays B — > 7r7r, contributions to the total amplitude from 
Ptttt and other amplitudes are comparable, such that the 
whole amplitude receives O(20%) corrections. This leads 
to corrections to the branching ratios and CP asymme- 
tries in B — ► 7T7T of order 



ABr(£ 
AA CP (B 



O(40%) 
O(20%) 



(58) 



These large uncertainties can be avoided by relying on 
isospin to define most of the parameters in the fit, as was 
done in Rcf. 

El 

in the e = method for 7 which has 
significantly smaller theoretical uncertainties. 

For B — > Ktt decays, the CKM factors and sizes of 
Wilson coefficients give an enhancement of the amplitude 
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parameter Pktt relative to the other amplitude parame- 
ters by a factor of order 10. Thus, the corrections to the 
total decay rates are suppressed by a factor of 10, while 
corrections to CP asymmetries, which require an inter- 
ference between Pkk with other amplitudes, remain the 
same. This gives 

ABr(B^Kn) ~ 0(5%) 
AA CP (B -» Kir) ~ 0(20%). (59) 

One exception is the CP asymmetry in B — > K°w~, 
which is strongly suppressed due to the smallness of the 
parameter A^. At subleading order A^-n can receive 
corrections far in excess of the leading order value, such 
that any numerical value of this CP asymmetry is com- 
pletely unreliable at the order we are working. We in- 
clude these estimates of power corrections into all our 
discussions below. 



IV. IMPLICATIONS OF SCET 

There are several simple observations one can make 
from the LO SCET expressions of the amplitude param- 
eters 

1. For C BMl ~ Cj Mi one finds that B -> MiM 2 
decays naturally have CmxM 2 ~ ^MiM 2 H2> so 
there is no color suppression. If one instead takes 
QBMt > QBMi as in Refs [gTJ then the u color 

suppressed" amplitude is indeed suppressed. 

2. There is no relative phase between the amplitudes 
CW, T^, T Kn , C K *, EW T K „ and EW C K „ and the 
sign and magnitude of these amplitudes can be pre- 
dicted with SCET. This allows the uncertainty in 
the Kir sum-rules to be determined, as well as pre- 
dictions for the relative signs of CP-asymmetries. 

3. The contributions of electroweak penguins, CV_io, 
can be computed without introducing additional 
hadronic parameters as discussed in Sec. IIII Dl 

4. The amplitude Aktc is suppressed either by A/m^, 
by small coefficients C^^, or by a s {m,b) compared 
with the larger Tktt and Cktt amplitudes 

5. If one treats (3k and /3 W as known, the amplitudes 

- c 

Tktt and EW K7T are determined entirely through 
the hadronic parameters describing the B — > itir 
system, implying that the branching ratios and CP- 
asymmetries for B — ► K + ir~ and B~ — > K°tt~ 
only involve 2 new parameters beyond 7T7t. 

6. In the combined SCET + SU(3) limit discussed in 
section IlII Bl the parameters — Pk-k — PrR? 



so we expect similar complex penguin amplitudes 
in B -> ATtt, B -> tttt, and B -> KK. 2 

Using these observations allows us to make important 
predictions for the observables, with and without per- 
forming fits to the data. Some of these have already 
been discussed and we elaborate on the remaining ones 
below. 



A. The ratio C/T and EW C /EW T 

We first describe in more detail the first point in the 
above list. Most literature has assumed that there is a 
hierarchy between the two amplitude parameters Cm 1 m 2 
and Tm 1 m 2 , ie. that Cm x m 2 ^ ^MiM 2 - This assumption 
is based on the fact that in naive factorization (in which 
C,/ = 0) one has C Mi m 2 /T Mi m 2 ~ c 2u /c lu ~ 0.1. The 
smallness of the ratio C2«/ci„ is due to the fact that the 
dominant Wilson coefficient C\ of the electroweak Hamil- 
tonian is multiplied by a factor of 1/N C in C2 U , explaining 
the name "color suppressed" amplitude, plus additional 
accidental cancellations which reduce the value of this 
ratio below 1/3. 

In SCET, however, the Wilson coefficients b\ : 2 con- 
tribute with equal strength to the overall physical ampli- 
tude and can spoil the color suppression 32]. In the bi 
terms for Cm-i_m 2 a factor of 1/N C = 1/3 occurs, however 
the hadronic parameter in the numerator is the inverse 
moment of a light cone distribution function and is ~ 3. 
Thus numerically Pn,K — 1, and setting (3 n> K = 1 for 
illustration we find 



T — 

± 7TTT 


-0.131C S7r 


- 0.099 Cj w 


Ctitt — 


-0.017^ 


- 0.160 Cj 7 " 


Tktt = 


-0.160 ( B7r 


- 0.120^ 


Crtt — 


-0.003 C S7r 


- 0.001 Cj" 




-0.0U( BK 


-0A59(f K 



(60) 



Thus, if £ ~ Cj it is easy to see that their is no "color 
suppression". On the other hand if ( 3> (j as chosen 
in Refs. [2^, |3(j then one would have significant color 
suppression. 

From Eqs. (|54I57(I the size of the color-suppressed and 
color allowed electroweak penguin amplitudes in 7T7t and 
Kir are directly related to that of CmiM 2 and Tm 1 m 2 - 
Thus if CmiM 2 ~ Tm 1 m 2 then SCET predicts that 
EWg lMa ~ EW£ lM3 . 



2 The analysis of "chirally enhanced" power corrections in Ref. |2-~f 
indicates that they will not break the equality in the SU(3) limit. 
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Parameter 


Measured value 


m B 


(5279.4 ± 0.5) MeV [60] 




(1.528 ± 0.009) ps [2] 




(1.643 ±0.010) ps [2] 


P 


0.379 ± 0.022 [2] 


u 


(130.7 ±0.4) MeV [60] 


Ik 


(159.8 ± 1.5) MeV [60] 


\Vud\ 


0.9739 ± 0.0003 [gjj 


\Vus\ 


0.2248 ± 0.0016 [61] 


\V a a\ 


0.2261 ±0.0010 [10] 


\Vcs\ 


0.9732 ± 0.0002 [10J 


\Vcb\ 


(41.6 ±0.5) x 10" 3 [2, 62] 


|Kb| incI 

1 t T 1 excl 
\Vub\ 

\V , |S IobaI 

|K 6 | hore 


(4.39 ±0.34) x 10" 3 [2, 63] 
(3.92 ±0.52) x 10" 3 [58, 64, 65] 
(3.53 ±0.22) x 10 -3 [10] 
(4.25 ±0.34) x 10 -3 



TABLE III: Summary of well measured input parameters. For 
our central value for \V u b\ we use a weighted average of the 
inclusive and exclusive |5S| with a slightly inflated error. 
Use ra t = 174.3 GeV. 



B. B — > 7T7T with Isospin and Im(C^ 7r /T7T7r) 







Using only SU(2) there are a total of 5 hadronic pa- 
rameters describing the decays B — > 7T7t, in addition to 
a weak phase. The 6 measurements allow in principle to 
determine all of these parameters as was first advocated 
by Gronau and London • Unfortunately, the large un- 
certainties in the direct CP asymmetry of B° — ► tt q tt° do 
not allow for a definitive analysis at the present time (ie. 
it currently gives 65° < a < 200° [13). It was shown in 
Ref. 0| that one can use SCET to eliminate one of the 5 
hadronic SU(2) parameters, since e = Im^C^^/T^^) ~ 0, 
and then directly fit for the remaining four hadronic pa- 
rameters and the weak angle 7, which substantially re- 
duces the uncertainty. Using the most recent data shown 
in section [J we find 



7" 



83.0°i^° ± 2 C 



(61) 



where the first error is from the experimental uncertain- 
ties, while the second uncertainty is an estimate of the 
theoretical uncertainties from the expansions in SCET, 
estimated by varying e ~ ±0.2 as explained in Ref. [4o| . 
This value is in disagreement with the results from a 
global fit to the unitarity triangle 



CKMfitter 
7 global 



58.6 



0+6.8 
5.9° i 



7ffi = 57 - 9 ° ± 7 - 4 ° 



(62) 



at the 2-ct level. A more sophisticated statistical analysis 
can be found in Ref. ^(|- The errors in Eq. I|61|) are 
slightly misleading because they do not remain Gaussian 
for larger e. At e = 0.3 the deviation drops to 1.5-er, 
and at e = 0.4 it drops to 0.5-cr. The result in Eq. I|61|) 



120 r 



Inclusive V u b 




direct 



3.25 3.50 3.75 4.00 4.25 4.50 4.75 



FIG. 1: Comparison of constraints on V u t and 7 from i) the 
direct measurement of /3, ii) current HFAG value for inclusive 
\V u b\, iii) global fit value of 7, iv) |Vut,| as output from the 
global fit |2lll0|. and v) results for 7 from the small e analysis 
of B — > 7T7T decays |45|. All errors bands are 1-a. 



is consistent with the direct measurement of this angle 
which has larger errors |2| 



7 



DK 



63 



o + 15° 
-12° 



(63) 



It is interesting to note that the global fit for (3 plus the 
inclusive determination of \V u b\ in table iLUl also prefers 
larger values of 7 as shown in Fig. ^ It will be quite 
interesting to see how these hints of discrepancies are 
sharpened or clarified in the future. In the remainder of 
this paper, we will show results for 7 = 83° and 7 = 59° 
to give the reader an indication of the 7 dependence of 
our results. 

The phase of the amplitude A™ is mostly determined 
from the CP asymmetries in B — * tt + it~. In particu- 
lar, as can be seen from the general parameterization of 
the amplitudes in Eq. (|13fl . the sign of the direct CP 



asymmetry C(ir^ 



is correlated with the relative sign 



between P n7T and and the sign of the asymmetry 
C(7T 7r°) with that between GV*- and P^-x- Since there is 
no relative phase between the amplitudes CV^ and T W7r 
at LO in SCET, the sign of the direct CP asymmetry 
in B — > 7r°7r° is thus expected to be positive based on 
the negative experimental value for C(n + n~) 46]. This 
expectation is in disagreement with the direct measure- 
ment shown in Table [I] Using the values of the hadronic 
parameters from the previous fit we find for 7 = 83° 



C(7r°7r°) = 0.49 ± 0.12 ± 0.23 , 
while for 7 = 59° we find 

C(7r°7r°) = 0.61 ±0.19 ±0.19. 



(64) 



(65) 



These values are 1.7cr from the measured value, if we add 
the theoretical and experimental errors in quadrature. 
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C. The decays B -> iyit in SCET 

For 7 = 83° a fit of the four SCET parameters to the 
B — > 7T7T data excluding the direct CP asymmetry in 
B — > 7r°7r° gives 



4.25 x 10~ 3 

|Kb| 
4.25 x ID" 3 



( Brr = (0.088 ± 0.019 ± 0.045) 
Cf T = (0.085 ±0.016 ±0.031) ( 



1Kb 



i(151±8±6)° 



IQ A P™ = (5.5±0.8±1.3)e 
while for 7 = 59° we find 

( Bn = (0.093 ± 0.023 ± 0.035) 

(f n = (0.10 ±0.016 ±.022) 

>3 TjUTt 



(66) 



4.25 x 10" 
I Kb I 
4.25 x 10" 3 
IKbl 



10 J P™ = (2.6±0.9±0.8)e 



i(103±19±16) c 



(67) 



The first error is purely from the uncertainties in the 
experimental data, while the second error comes from 
adding our estimate of the theory uncertainties discussed 
in section IIII El For both values of 7 one finds that 
\P/T\ ~ 0.25. Note that this ratio of P/T does not 
include the ratio of the CKM factors. The ratio relevant 
for the decays is \>& /)$\\P/T\ ~ 0.6. 

It is interesting to compare the result for P 7 ™ extracted 
from the data with that from the purely perturbative 
penguin computed in Ref. [25j (in two scenarios for the 
input parameters), 



1 n3 r>dcfault 
iU -fQCDF 



-1.0 - [0.1 + 0.3« + . . .] (68) 

- (1.7 + Q.QXh + .047X1) - {0.3 + 0.H} 

= -3.1 -0.4?: -0.047^1, 
10 3 -Pqcdf = -0-9 - [0.0 + 0.1* + . . .] 

- (2.0 + 0.0X H + .063X1) - {0.3 + 0.M} 

= -3.2 - OAi - 0.063X1 , . 

Here the terms [. . .] are a s (mb) corrections, the terms 
(...) are chirally enhanced power corrections with param- 
eters Xh and Xa, and {. . .} are perturbative corrections 
to these. We observe that the magnitude of the pertur- 
bative 7T7T penguin is of similar size to that from the data 
for 7 = 59°, but has a small strong phase is in contrast 
to the large strong phase seen in the data. 

The correlation between C, B * and Cj w in Eqs. (|66I67|) is 
about —0.8, so that the heavy-to-light form factor, which 
is given by the sum of these two parameters is determined 
with much smaller uncertainties than one would obtain 
by naively adding the two individual errors in quadra- 
ture. For 7 = 83° we find 



F BlT {0) = (0.17 ± 0.01 ± 0.03) ( 



4.25 x 10" 3 



(69) 



while for 7 = 59° we find 

F B7! (0) = (0.19 ± 0.01 ± 0.03) ( 



4.25 x 10" 

IKJ 



(70) 



D. The decays B -> Kir 



For these decays the penguin amplitudes are enhanced 
by the ratio of CKM matrix elements \\ { c s) /\i s) \ ~ 40. 
Thus, the relevant ratio of penguin to tree amplitudes 
is |Ai s) /Al s) ||P/T| - 10 and the B -> Ktt decays are 
penguin dominated. If one were to only keep the pen- 
guin contributions to these decays the relative sizes of the 
branching ratios would be determined by simple Clebsch- 
Gordon coefficients 



Br(7r u in ~ Br(7r u ir 



(71) 



Deviations from this relation are determined at leading 
order in the power counting by the non-perturbative pa- 
rameters ( BM and Cj M ■ To see how well the current 
data constrains deviations from this result we can look 
at the following ratios of branching fractions 



Ri = 



R2 



2Br(£- -> 7r°K~-) 



Br(B- 
Br(B° 



-> TT-K°) 

tt-K+)t b 



- 1 = 0.004 ± 0.086. 



(72) 



R3 



Bt(B- 
2Bt(B~ 



7T K°)tbO 

, 7r°K )r B - 



- 1 = -0.157 ±0.055. 



1 = 0.026 ±0.105. 



Br(B° -> tt-K )t b o 

and the rescaled asymmetries 

Ai = (1 + Ri)A C p(n°K-) = 0.040 ± 0.040 , (73) 

A 2 = (1 + R 2 )A C p(n-K + ) = -0.097 ±0.016, 

A 3 = (1 + i? 3 )A C p(7r X ) = -0.021 ± 0.133 , 

A 4 = ^cp(7t"X ) = -0.02 ±0.04. (74) 

These ratios have been defined by normalizing each 
branching ratio to the decay B~ — * ir°K°. If we drop 
the small amplitude parameter Ak-k then this channel 
measures the penguin, 



A(B~ ^n-K°) = X^P Kw 



(75) 



and the direct CP asymmetry is expected to be small. 

A simple test for the consistency of the Kw data is 
given by the Lipkin sum-rule for branc hing ratios |55| . 
and a sum-rule for the CP-asymmetries 69] 



Ri - R2 + Rs = , Ai 



0, (76) 



which are both second order in the ratio of small to large 
amplitudes as discussed in section lTl El The current data 
gives 

R 1 -R 2 +R 3 = (0.19 ± 0.15) oxpt , 

A x - A 2 + A 3 - A 4 = (0.14 ± 0.15) cxpt . (77) 
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Thus, so far this global test does not show a deviation 
from the expectation. 

SCET provides us with additional tests for the Kir 
data. It turns out that the current data is not precise 
enough to determine the values of £ BK and ( B K . These 
two parameters only contribute to the two decays B~ — > 
%°K~ and B° — > ir°K°, which have neutral pions and 
larger experimental uncertainties. As we will explain, 
the data on these decays seems to favor a negative value 
of C,j K , but that would imply a negative value for X B , 
the first inverse moment of the B meson wave function, 
contrary to any theoretical prejudice. One can use the 
fact that the only sizeable strong phase is in the value 
of the parameter A^J to determine the predicted size of 
the deviations from the above relations and also the signs 
and hierarchy for the CP asymmetries. 



1. Sum-Rules in B — > Kn 



In SCET positive values of ( BM and (f M imply that 
the phase of -T Kw , -C Kw , EW G £, and EW% - EW%„ 
are the same. This can be seen from Eas. (l49l5UII . There- 
fore this implies that these amplitudes have a common 
strong phase S relative to the penguin Pk-k- Using the 
notation and results from section IlIEJwe have 



we find 



ew 



1.40(C S7r + C,D + 0.35%Cj', 

0.12(C bjc +Cj K ) + 1.27AtCj*, 
+0.03(C S7r + C.f *) - Q.02#jcCj* 
0.7l(( BK + (f K ) - O.170„{f K 

+0.12(C S " + C.f 7r ) + 0.87/foCj w 

0.12(C S " + C,f") + 0.87ArCj* • 



Generically Q BM + Cj M 



0.15 



0.25 and (f M 



(84) 
0.05- 



0.15 so that 6t, ec, e^J", eg" are ~ 0.1-0.4 and can be 
thought of as expansion parameters. 

To estimate the SM deviations from the results in 
Eq. O we take the C(e 2 ) terms in Eqs. (|30I31[) and 
independently vary the parameters in the conservative 



ranges ( B7T + ( B 71 



0.2 ± 0.1, f3 R ( B n = 0.10 ± 0.05, 



0.2±0.1, ArCf* 



s , - 0.10±0.Q5, e A = 0±0.1, 

7 = 70° ±15°, arbitrary <^4 and all phase differences 
A^> = 0° ± 30° . For the Lipkin sum rule this gives 



R 1 -R 2 + R 3 = 0.028 ± 0.021 , 
and for the CP-sum rule 

Ax - A 2 + A 3 - A 4 = ± 0.013 . 



(85) 



(86) 



be = <t> : 



6f 



(78) 



At LO in SCET one can drop the Ak-k amplitude (e^ 
0) and write 



A(B° - 




A(B° - 


> 7T+K-) 


V2A(B- - 


■+ TT°K-) 


V2A(B° - 





(79) 



l + T (^-^e-^) 



~2~ 



where the e-parameters are all positive and satisfy 



. , j_ , ,ew ew ew (or\\ 

e — e T + e c , e — e T — e c , (8U) 



Experimental deviations that are larger than these would 
be a signal for new physics. The CP-sum rule has sig- 
nificantly smaller uncertainty than the Lipkin sum-rule. 
This can be understood from the expression 



Ai - A 2 + A 3 - A 4 = - 

— £t^c" sin((f>T 
x (l + 0(e A )) . 



sin( 7 ) [4™ esin(<£ 



e c e ew sm 



(4> C - 4>ew)] 



(87) 



All terms involve one of the smaller electroweak penguin 
e-parameters, and in SCET all the phase differences are 
small, both of which give a further suppression over the 
Lipkin sum-rule. Since the CP sum-rule is always sup- 
pressed by at least three small parameters it is likely to 
be very accurate. 



and 



e > e c , e > e T . 



(81) 



From the decomposition in terms of SCET parameters 
we can determine the magnitudes of the e-parameters in 
terms of the £'s. The rate Br(_B _ — > ir~K°) determines 



W 3 \P K J ~ 5.5 ±0.1 ±0.1 



(82) 



2. CP-Asymmetry Sign Correlations 

For the asymmetry parameters up to smaller terms of 
0(e 2 ) we have 



Ai = — esin((5) sin(7), 
A 2 = — er sin((5) sin(7), 
A 3 = e c sin(<5) sin( 7 ). 



(88) 



and using 



Thus, we immediately have the following predictions 



0.0236, 



(83) 



i) Ai, A 2 , —A 3 have the same sign , 

ii) |Ai|>|A 3 |, |A 1 |>|A 3 |, (89) 
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where i) depends only on the fact that positive £'s 
gives positive e-parameters, and ii) follows from including 
Eq. ijHH - Compared to the data in Eq. if?^ we see that 
the central values of Ai and A2 currently have opposite 
signs, disagreeing from equality by ~ 2a when we take 
into account the theoretical uncertainty. The experimen- 
tal errors are still too large to draw strong conclusions. 

Note that a prediction JAi| w IA2I was made for the 
CP-asymmetries in Ref. [5(j based on the expectation 
that the color suppressed amplitudes are small. The CP- 
sum rule Ax — A 2 + A3 = was discussed in Ref. j6^ 
(3rd reference) to take into account the possibly large 



color suppressed contributions. Given £ 



Cj , SCET 



predicts that the phase of the color suppressed Ck-k am- 
plitude is nearly equal to that of the Tk-k amplitude so 
the hierarchy of the asymmetries is actually reinforced 
by a significant Ck-k- Our prediction that |Ai| > |A 2 | 
with Ai ; 2 having equal signs can also be compared to pre- 
diction for the analagous CP-asymmetries in the QCDF 
approach [2^] (4th reference). Four different scenarios 
for the hadronic parameters were considered S1,S2,S3,S4, 
and all four sets of model parameters exhibit the sign cor- 
relation. (However all four of the scenarios also underes- 
timate the size of |Acp(7r + if _ )| by more than a factor of 
two due mostly to the fact that the purely perturbative 
penguin for Kn is somewhat small.) 

For the branching ratio deviation parameters we have 
up to smaller terms of 0(e 2 ) that 



Rt = cos((5)[e^-ecos(7)], 
R 2 = cos( ( 5)[eg u -e T cos(7)], 
R 3 = cos(<5)[e c cos( 7 )-e e ™l. 



(90) 



The use of conservative errors on the ^-parameters leaves 
too much freedom to make sign predictions for the R^s. 
However, definite sign predictions will be possible us- 
ing Eq. H90|) when the £ parameters are pinned down 
by B — > n and B — > K form factor results in the fu- 
ture. Alternatively accurate measurements of the Ri 
plus Acp(K + tt~) will determine the hadronic parame- 
ters needed to predict the magnitude of the remaining 
Ai's. 

In the next section we turn to more direct comparisons 
of the SCET predictions with the data by fixing the pa- 
rameters with the well measured observables and then 
predicting the rest. 



B- 



n~K° and B -> ir~ iT 1 



The amplitude parameters Tk-k, Ak-k and EW Kn are 
determined in terms of the parameters C BTr and (f n ob- 
tained previously from the decays B — > nn. Thus, only 
two new parameters are required for the decays B~ — > 
n~K° and B Q — ► n + K~: the magnitude and phase of 
P K7r . Since the ratio of X u A K7r < \ c P K jr ~ 0.00L, one 
predicts a negligible CP asymmetry in B~~ — > n~ K° in 
agreement with the data. The best sensitivity on the two 



parameters is from Br(£> — > n K°) and Acp{B° — > 
n + K~). Using these two observables we find two solu- 
tions for A** for 7 = 83° 



(91) 



1 n 3 pK-K _ I (5.5±0.1±0.1)e*( 144±8±11 >° 
\ (5.5±0.1±0.1)e 4 ( 32±7±1 °)° 

while for 7 = 59° we find 

m SpK«_i (5.5±0.1±0.1)e^ 44±9±11 )° 

\ (5.5±0.1±0.1)e*( 36±8±10 )° 1 ' 

The confidence level plot for the magnitude and phase 
of P K7r is shown on the left of Fig. El For the 7 = 59° 
result the magnitude indicates a large SU(3) violating 
correction at leading order in A/E^ or a large A/E^ cor- 
rection in the SU(3) limit (which disfavors this solution). 
Taking the 7 = 83° we see that of the two solutions the 
first has a phase which agrees well with the SU(3) rela- 
tion to the phase in 7T7t, while the second phase is quite 
different. 

For 7 = 83° the first solution, however, does not give 
good agreement with the third piece of data, the branch- 
ing ratio Br(B° ->■ tt+K~) = (18.2 ± 0.8) x 10" 6 , while 
the second agrees considerably better. We find 



Br(7r+Jr 



(24.0 ±0.2 ±1.2) x 10~ 6 
(21.3 ±0.2 ±1.3) x 10~ 6 



(93) 



For 7 = 59° this branching ratio has much less discrimi- 
nating power between these two solutions and we find 



+ f (22.5 ±0.2 ±1.2) 

I (22.7 ±0.3 ±1.2) 



x 10~ 6 
x 10~ 6 



(94) 



This is can also be clearly seen in the confidence level 
plot for P K7r on the right of Fig. |2 where we have in- 
cluded the branching ratio measurement in the fit. Note, 
however that both solutions have trouble explaining the 
small branching ratio Br(_B° — > ir + K~), making the 
large difference in the branching ratios of B — > ir + K~ 
and B — > ir~ K° quite difficult to explain at LO in the 
A/m b < 1 limit of QCD. 



4- Predictions for other Kn and KK observables 



Using the hadronic parameters extracted from the 
B — > 7T7T decays (( B7r , (f^ and P„), the value for Pk^ 
determined from the decays B~ — > n~K and B° — > 
it~K + decays and independently varying C, BK + (f K = 
0.2 ± 0.1 and /3 w (f K = 0.10 ± 0.05, we can calculate all 
the remaining currently measured Kn observables. The 
results are given in Table Hvl for 7 = 83° and 7 = 59°, 
respectively. We also show these results in Figs. |3 and 0] 
The data used in the fit are shown in red below the 
dashed dividing line while those above the line are pre- 
dictions. Note that there is one more piece of data below 
the line than there are parameters. 
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Penguin in the Kn system 




4.5 4.75 5.0 5.25 5.5 5.75 6.0 6.25 4.5 4.75 5.0 5.25 5.5 5.75 6.0 6.25 



Penguin in the tie system 

y = 79 y =59 

\P^\ =5.5 ±1.5 \P^\=2.6±1.1 
arg(P^) = 151°+ 10° arg(P^) = 103°±25° 

FIG. 2: Confidence level plots for the complex parameter A^J 
for 7 = 83° (left hand side) and 7 = 59° (right hand side). On 
the top we show the confidence levels without using Br(B — > 
ir + K~), while the bottom plot includes this branching ratio. 
We also show the value of P,™-, which is idential to the Pktv 
in the SU(3) limit. 





Expt. 


Theory 


Theory 






(7 = 83°) 


(7 = 59°) 


Data in Fit 




sv+o 


-0.50 ±0.12 


-0.50 ±0.10 


-0.51 ±0.10 


C(7r + 7r-) 


-0.37 ±0.10 


-0.37 ±0.07 


-0.38 ± 0.07 


Br(7r + 7T~) 


5.0 ±0.4 


5.0 ±2.0 


4.6 ± 1.8 


Br(7r+7r°) 


5.5 ±0.6 


5.5 ±2.2 


7.3 ±2.9 


Br(7r°7r°) 


1.45 ±0.29 


1.45 ±0.58 


1.32 ±0.53 


Br(^ 7O 


24.1 ± 1.3 


24.1 ± 1.2 


24.1 ± 1.2 


A{K-n+) 


-0.115 ±0.018 


-0.115 ± 0.023 


-0.115 ±0.023 


Br(K°K-) 


1.2 ±0.3 


1.2 ±0.5 


1.2 ±0.5 


Predictions 




A{tV + 7T°) 


0.01 ±0.06 


< 0.05 


< 0.05 


A(7V°7T°) 


0.28 ± 0.40 


-0.48 ±0.19 


-0.52 ±0.27 






0.84 ±0.23 


-0.14 ±0.22 


Br(ir°K°) 


11.5 ± 1.0 


10.4 ± 1.1 


10.9 ± 1.2 


Br( 7 r+A'-) 


18.9 ±0.7 


24.0 ± 2.1 


22.5 ± 2.1 


Bt(it°K-) 


12.1 ±0.8 


14.3 ± 1.5 


12.7 ± 1.4 


S(rr K s ) 


0.31 ±0.26 


0.77 ±0.16 


0.76 ±0.16 


A(ir°K-) 


0.04 ± 0.04 


-0.183 ±0.075 


-0.184 ±0.076 


A(K°w°) 


-0.02 ±0.13 


0.103 ±0.058 


0.083 ± 0.047 


A(tt-K°) 


-0.02 ±0.04 


< 0.1 


< 0.1 


Br(K°K°) 


0.96 ±0.25 


1.1 ±0.3 


1.1 ±0.3 


Br (K + K~) 


0.06 ±0.12 


<o.i 


< 0.1 


A(K°K~) 




<0.2 


^ 0.2 


A(K°K°) 




<0.2 


^ 0.2 



TABLE IV: Comparison of LO predictions versus data as in 
Figs. 3,4. Br's are in units of 10~ 6 . The theory errors dis- 
played for quantities used in the fit show the relative weight 
for these observables from power corrections, while those for 
predictions include parameter uncertainty from the fit as well 
as from power corrections. CP asymmetries that are not 
shown in the table are not determined at this order. 



In Fig. we see that 7 = 83° gives a good match to the 
B — > 7T7T data except for the asymmetry C(7r°7r°). When 
taking into account the theoretical error the most strik- 
ing disagreements are the Br(K~ir + ) at 2.3er and the 
CP-asymmetry Aqp(K ~7t°) at 2.6a. All other predic- 
tions agree within the uncertainties. Note that one could 
demand that Acp(K~ir°) be reproduced, which would 
imply a negative value of £j K (a naive fit for 7 = 83° 
gives (,f K ~ —0.15). Note however, that this would im- 
ply that both perturbation theory at the intermediate 
scale fi = VEA and SU(3) are badly broken. 

The situation in Fig. 0] with 7 = 59° is similar except 
that the theoretical prediction for Br(7r + 7r°) moves some- 
what. The Br(K~ir + ) deviation is reduced to 1.6a and 
asymmetry A(K~tt°) is still 2.6a. All other predictions 
agree within the uncertainties. 

For B — > KK the amplitude parameters in SCET sat- 
isfy A KK = B KK and E KK = PA KK = EW KK = 0, 



and we obtain the prediction 

Br(fl _ -> K-K°) = Bt{B° -> K°K°) (95) 

which agrees well with the latest data, and the expecta- 
tion that A CP (B~ -> K~K°), A CP {B° -> K°K°), and 
Br(£>° — > K~K + ) will be suppressed. 

Unfortunately, without the use of SU(3) we do not 
have enough experimental information to determine the 
hadronic parameters required to predict the B — > K°K° 
absolute branching ratio. It is however interesting to ex- 
tract the penguin amplitude and compare with the other 
channels. Wc find 

W 3 \Pkk\ = 5.3 ±0.8 (96) 

Comparing with the penguin amplitudes extracted in ttt: 
and in Kn we see that the combined SU(3) and SCET 
prediction, P^^ ~ P Kt . ^ Pkki works quite well if 7 = 
83°. 
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The CP asymmetries 



The CP asymmetries 



S(7l°7l°) 
A(7t°7C°) 

A(K 71°) 
A(K° %-) 
A(K°7l°) 
S(K°7l°) 
A(K"7i+) 

C(7l+7l-) 
S(7t + 7T) 



K°K° 

KV 
K"7l° 
K 7C + 

K tT 
K°K° 

71° 71° 
7C + 7C° 
7t + 7lT 



. 1 A 1 , 


i * 1 




l * 

■ .— #— i 


* 1 

* y=83° 

1 * ! 1 4 ' 




A Theory 

* Data 

. ... i .... i .... i .... i . . 



-0.5 -0.25 0.25 0.5 0.75 i.O 
77te Branching ratios (xlO 6 ) 




A 



•—In 

I 

5 



Theory 
Data 



S(7C°7C° 
A(7t°7t° 
A(K 7C° 

A(K° 7T 
A(K°7t () 

S(K ( V 
A(K"7C + 

C(7C + tc- 

S(7t + 7T 




The Branching ratios (xlO ) 



10 



15 



20 



25 




7C + 7C 



FIG. 3: Comparison of theory and experiment for all available 
data in B — > tttt and B — > i^7r decays, with 7 = 83°. The 8 
pieces of data in red (below the dashed line) have been used 
to determine the SCET hadronic parameters C, BlT , Cj^j Pmr, 
Pk-tt and \PkK\, with ^" s/f and £f K fixed as described in 
the text. The data above the line are predictions. The CP 
asymmetry in B~ — > K°ir~ is expected to be small, but its 
numerical value is not predicted reliably. 



V. CONCLUSIONS 

Decays of B mesons to two pseudoscalar mesons pro- 
vide a rich environment to test our understanding of the 
standard model and to look for physics beyond the stan- 
dard model. The underlying electroweak physics mediat- 
ing these decays are contained in the Wilson coefficients 
of the electroweak Hamiltonian as well as CKM matrix 
elements. In order to test cleanly the standard model 
predictions for these short distance parameters, one re- 
quires a good understanding of the QCD matrix elements 
of the effective operators, which can not be calculated 
perturbatively. 

At the present time, there are 5 well measured (with 
< 100% uncertainty) observables in B — > tttt, 5 in 
B — > Ktt and 2 in B — > KK . Using only isospin symme- 
try (with corrections suppressed by m Ui d/A), the number 
of hadronic parameters required to describe these decays 
is 7, 11 and 11, respectively. The number of hadronic 
parameters can be reduced by two in the rnr system, 



FIG. 4: Same as Fig.3, but with 7 = 59°. 



if one drops the two operators O7 and O s , which have 
small Wilson coefficients in the standard model. If one 
is willing to take SU(3) (an expansion in m s /A) as a 
good symmetry of QCD, the combined B — > tttt /Ktt sys- 
tem is described by 15 parameters, while the B — > KK 
system adds another 4 parameters. Neglecting O7 and 
Os with SU(3) reduces the number of parameters in the 
tttt I 'Ktt I KK system to 15. Thus, at the present time 
there are more hadronic parameters than there are well 
measured observables. 

In this paper we have studied these decays in a model 
independent way using SCET. This analysis exploits that 
the hadronic scale A in QCD is much smaller than both 
in the large mass of the heavy quark and the large energy 
of the two light mesons. It follows that at leading order 
in the power expansion in Aqcd/Q, where Q ~ rrib,E, 
and using SU(2), there are four hadronic parameters de- 
scribing B — > 7T7T, five additional parameters describing 
B — > Ktt and three additional parameters describing 
B — * KK. In the limit of exact SU(3) the four param- 
eters describing B — > tttt are enough to describe all of 
these B -> PP decays in SCET. 

In SCET the electroweak penguin operators O-j^ can 
be included without adding additional hadronic param- 
eters. One can use the 5 pieces of well measured 7T7T 
data to determine the four hadronic parameters and the 
weak angle 7 |45j . and with the current data one finds 
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7 = 83° ± 8° ± 2°. This is still consistent with the direct 
measurement of this angle from B — ► DK [2| , but is cur- 
rently in conflict with the value of 7 from a global fit of 
the unitarity triangle at the 2a level. It is too early to tell 
if this implies larger than expected power corrections in 
SCET or might be a first hint at new physics. When we 
proceed to analyze the decays B — > Kn, we thus perform 
our analysis both for 7 = 83° and 7 = 59°. For both of 
these values the direct CP asymmetry in B — > 7r°7r° is 
predicted to have the opposite sign from the measured 
value, but is still consistent at the 2a level. 

Moving on to B — ► Kir decays, we analyzed the uncer- 
tainty in the Lipkin sum rule |55j for branching fractions 
and the CP-sum rule j6(| for rescaled CP-asymmetries as 
defined in Eqs. H72I73|) . giving our result in Eqs. (|85I86|) . 
The CP-sum rule was found to be particularly accurate 
due to a suppression by three small parameters in SCET. 
The Lipkin sum rule is second order in small parameters 
and has a theoretical precision that also makes it an in- 
teresting observable. We conclude that both the Lipkin 
and CP-sum rules will provide very robust methods for 
testing the Kit data as the experimental errors decrease 
in the future. 

Using the expectation that the hadronic parameters 
<^ BAI and Cj m the factorization theorem are posi- 
tive, we showed that the rescaled asymmetry Ai(ir°K~) 
should have the same sign and larger magnitude than the 
rescaled asymmetry A2(7r~ K + ) which is well measured. 
This prediction is in conflict with the current data by 
~ 2cr. Other sign and magnitude predictions are dis- 
cussed in section flV D 21 

The SCET amplitude formulas predict that in addition 
to the 7T7T parameters already determined, only the com- 
plex Kit penguin amplitude is required to describe the 
decays B~ — ► ir~K° and B — ► ir + K~. This happens be- 
cause they involve ( B7V and Cj ir , but do not involve C, BK 
or Cj K ■ The well know prediction of a small CP asymme- 
try for B~ — ► tt~K° is reproduced in SCET. The large 
difference in Br(B — ► ir + K~) and Br(i?~ — > ir~ K°) is 
difficult to explain in the standard model with SCET. 
The 7 = 59° solution is not preferred by the combined 
SU(3)+SCET limit which predicts P Kn ~ P^. These 
amplitudes agree well for 7 = 83°. 

Given the current uncertainties in the data, the re- 
maining two hadronic parameters C, BK and Cj can not 
yet be determined reliably. This also means that pre- 
dictions for the remaining rates do not depend too sensi- 
tively on these parameters. Fixing their values to be close 
to those preferred by SU(3), but with 50% uncertainty, 
we obtained predictions for the remaining observables in 
Figures and 0] 

Finally, the decays B — > KK require two additional 
hadronic parameters, which can only be determined once 
better data for both rates and CP asymmetries become 
available for these decays. One prediction of SCET, 
namely that Br(S -> K°K°) = Br(B~ -> K~K°) is 
well satisfied by the current data. In the SU(3) limit 
one expects that Pk-k ~ PrRi an d this result is in good 



agreement with the data. 

In conclusion, several predictions of SCET work rather 
well, while for others there are discrepancies with the 
current data. It is too early to tell if the disagreements 
between theory and data are due to statistical fluctu- 
ations, to larger than expected power corrections or if 
they reveal a first glimpse of physics beyond the stan- 
dard model. To answer this question, the experimen- 
tal uncertainties need to be reduced and the convergence 
of the SCET expansion of QCD for nonleptonic decays 
has to be tested further both with nonleptonic and with 
semileptonic data |65ll70j. 
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APPENDIX A: OPERATORS AND MATRIX 
ELEMENTS IN SCET 

At the scale /x ~ m,b the Hamiltonian in Eq. JSJ is 
matched onto operators in SCET. For the first two orders 
in the power expansion 

* n.n ^ i 

[duAUWMQW&i) + + ■ ■ • } . ( A1 ) 

The Wilson coefficients Cj and 6, are the Wilson coeffi- 
cients that appear in Eqs. I|34II36II . The operators for the 
AS = transitions are |2|| H3] 

Qld = [Wn,uu [dn,u 2 iPLUn^ 3 ] , (A2) 
Qfdfid = [dn,Ui ^PlK] [q.n,u 2 i-PL,Rqn,u 3 \ , 
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and 



TTl}) 

2 

TTl}) 

2 

mi, 
2 

i 'lb 



o (1) 

^7d 



RUf, 



Tflb 

The AS" = 1 operators Q\g are obtained by swapping 
d — > s. The "quark" fields with subscripts n and h are 
products of collinear quark fields and Wilson lines with 
large momenta u>{. We have defined 



1 



2 9 



(-«) 



(A3) 

[Wl[in-D Cin , %D» ± }W n 6(uj - pt)] 



where creates a collinear up-quark moving along the 
n direction, or annihilates an antiquark. The b v field is 
the standard HQET field. For a complete basis we also 
need operators with octet bilincars, T A <g> T A , but their 
matrix elements vanish at LO. The operators Q^J and 

Qg) also do not contribute at LO [12, see also |39| . 

The leading order factorization theorem in Eq. I|32() is 
generated by time ordered products of both the operators 
and with insertions of a subleading Lagrangian. 
T-products with contribute to terms with Q BM and 
T-products with Q^> contribute to those with (f M . It 
is convenient to define 



Q 



(i) 



nib 



Qt = qI^PlMI 



(A4) 



(A5) 



In Q^ ' 1 ^ the flavor of the q\ u terms matches that of the 
first bilinear in Eq. (|A2I) . In Q™ the flavor of q 1 and q n 
match those in the second bilinear of Eq. I|A2(1 . and we 
have Pr for i = 3, 6 and Pl otherwise. The contributions 
to B — > Mi Ma at LO are all from Q™ times the time- 
ordered products 

T{=jd±yT[QP (0) 

+ JdVV r[Qf (o) *4!>) 
+ /Ar[g l (0) (o), ? 4 1 „f( 2 /)], 

=/AT[Q l (1) (0),z4 1 j g ( y )], (A6) 



where z and 1 — z are the momentum fractions carried 
by the collinear quark and gluon field in Q^p ■ Here 
T\ and are exactly the same T-products that occur 
in the heavy-to- light form factors [66j. In addition we 
have operators/T-products whose matrix elements give 
A cc (see the appendix of Ref. [34[ for further discussion 
of these contributions). Using the collinear gluon fields 
defined in Ref. the Lagrangians in Eq. (|A6|) are 



4? = (^W) l ^ s i(wUpt^n)+h.c, (A7) 



- £ — ^ 



(2) 



$ 1 



igftyWqus + h.c. 



2 m-D, 

+Cn I 777-7777 «S^VFg us + h.c, 



2 ^{in-Dc) 2 
4^ = ^tr{[iD$,iD±»][iD 0lt) WiDt. v Wi]} , 



where i£>£ = ilX 1 + gA£. 

In this paper we only used this factorization at the 
scale rrib, so the hadronic paramaters are defined by ma- 
trix elements of T\ and Ti and the n-collinear operator, 
namely 

(M n \Tt\B) = d(B, M) m B Q BM , (A8) 
(M n \T*(z)\B) = d(B, M) m B (f M (z) , 
(Mn\Qf |0) = Cl(B, M) m B fa cj) M (u) 

where u and 1 — u are momentum fractions for the 
quark and antiquark n-collinear fields. Here C l (B,M) 
and C[(B,M) are simple Clebsch-Gordan coefficients. 
Putting the pieces together we have 

A = (M 1 M 2 \H W \B) (A9) 

= J2 C ^ B > C '^ B ' M 2) fM 2 



dudz bi(u, z)Cj Ml (z)4>m 2 ( u ) 



+ C BMl / du Cl {u)$M 2 (u) 
Jo 



+ (!«-> 2) 



This result was used to obtain Eq. (|34H3fci|) where the 
relevant combinations of C l C n coefficients can be read 
off from Table I of Ref. [13 (and do not asssume isospin 
symmetry). Here A c ^ 2 contains Clebsch Gordan coef- 
ficients if for example SU(2) is used to relate these pa- 
rameters in different channels. For amplitudes with no 
penguin contribution we have A^ M2 = 0. 
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APPENDIX B: RELATIONSHIP BETWEEN OUR 
AMPLITUDE PARAMETERIZATION AND 
GRAPHICAL AMPLITUDES 

In this appendix we show the relationship between the 
amplitude parameters defined in Eqs. JT3J-|JTSJ| and the 
graphical amplitudes defined in |23Ll47| . These relations 
are useful, since one can immediately read off SU(3) re- 
lations between different amplitudes, since the graphical 
amplitudes are SU(3) invariant. Note that while the am- 
plitude paramters on the right hand side of equations 
HB1IB21 153)1 have the same name for the different pro- 
cesses, 7T7T, Kit, and KK, they are only equal in the 
SU(3) limit. 

The relations for the amplitude parameters in B — > txtx 
are 



= T + P ut + E + PA ut + EW C 
EW P EW PA 



EW EW 



E 



C-Tr-n- — C — P u t — E — PA u t 



3EW T EW C 



2 2 
EW P EW PA EW E EW A 



2 2 2 

EW A EW E 



written in terms of graphical amplitudes as follows: 

EW P EW E 



P™ = Pct + PA ct + EW' 

2 A 

EW P EW PA 
2 2 

EW1„ = ^(EW T + EW C ). (Bl) 
The amplitude parameters for B — > Kit decays can be 



T Kv = T + P ut + EW c - 



Cktt = C — P, 



2 2 
3EW T EW C 



2 2 
EW P EW E 



Pk* = P ct + EW E - 



EW C EW P 



A K n = Put + A + EW - 

EW T K7T = ^(EW T + EW C ) 

EW° Kn = ^(EW C -EW E ), 
Finallyfor B — > KK decays we find 



2 2 

EW P EW C 



(B2) 



A K k = Put+A- 



EW 



c 



EW - 



EW 1 



B KK = Put + PAut - 



Ekk — —E — PAut — 



EW 



- EW J 



EW A EW PA 



Pkk = Pct- 



EW 



c 



EW - 



EW 1 



PA KK = PA, 



2 2 
EW A EW PA 



EW 



2 2 
3EW A 3EW E 



KK — 



(B3) 
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